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ABSTRACT 

The LINAC at the Naval Postgraduate School , Monterey, was used 
to accelerate electrons to energies ranging from 52 to 92 MeV in 
order to study the energy distributions of high energy electrons 
before and after passing through layers of tin, gadolinium, and 

2 

lead. The thickness of these materials ranged from 0.8 to 5.9 g/cm . 

The most probable energy losses agreed with the theory of Blunck 
and Westphal for all materials used, while distribution half-widths 
agreed only for absorbers of thickness less than 3.0 g/cm . The 
thickness at which theory and experiment began to exhibit a notice- 
able discrepancy was found to be dependent on the atomic number of 
the material. 

Where comparison was possible, results of this experiment generally 
agreed with the findings of similar works concluded previously. 
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I. INTRODUCTION 



The theoretical energy distribution for an initially monoenergetic 
electron beam which has passed through an absorbing material was de- 
termined by Blunck and Westphal [1]. Their work modified previous 
theoretical calculations of Landau [2], Eyges [3], Bethe and Heitler 
[4], and Blunck and Leisegang [5] by adding radiation losses to the 
previous predictions of ionization losses for electrons when passing 
through thin absorbing layers. The distribution assumes energy losses 
small compared to the incident beam energy. The theoretical treatment 
is presented in Section II. 

Several measurements of energy losses by high energy electrons 
passing through various materials have been performed in the 10 MeV 
to 150 MeV range. Of note are the works of Breuer on aluminium [6]; 
Bumiller, Buskirk, Dyer, and Miller on aluminum [7,8]; Goodwin on 
copper [9]; and DeLeuil and Raynis on aluminum, copper, and lead [10]. 
Except for Breuer' s work on the Darmstadt linear accelerator, the 
experiments were performed on the LINAC at the Naval Postgraduate School . 
These experimental results are in general agreement with theory for thin 
absorbers (s 2 g/cm ), but discrepancies with theoretical predictions 
are reported for the thicker absorbers. Such discrepancies are expected 
for the thicker absorbers since the energy loss is no longer small 
compared to the initial beam energy. 

In this thesis the energy loss measurements are extended to include 

tin (Z = 50) and gadolinium (Z = 64). These metals were chosen to 
provide experimental results tor materials with atomic numbers between 
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copper (Z = 29) and lead (Z = 82). It was thought that analysis of 
these materials would provide more conclusive information on the 
effects of atomic number on the agreement between experimental and 

theoretical half-widths as reported by DeLeuil and Raynis [10]. 

2 

Thicknesses used ranged from 1.485 to 5.940 g/cm for tin and from 
0.814 to 4.831 g/cm for gadolinium; nominal beam energies were 52, 

75, and 92 MeV. Thick lead absorbers (2.825 and 4.236 g/cm^) were 
investigated as an extension of the work of DeLeuil and Raynis, and, 
in addition, their thickest aluminum absorber (5.574 g/cm ) was 
investigated as a continuity check between the two experiments. 

The energy loss distributions are characterized by the most probable 
energy loss and the half-width. The half-width is the full width of 
the distribution curve at half maximum. These quantities, obtained 
from experiment and theory, are the basis for comparison between ex- 
perimental and theoretical results. 
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II. THEORETICAL CONSIDERATIONS 



The Blunck and Westphal theory of the distribution for the energy 
loss of a beam of monoenergetic electrons in passing through a layer 
of absorbing material assumes that the energy loss Q is small compared 
to the initial beam energy, E^. Let W(Q)dQ be the probability of 
energy loss between Q and Q + dQ, and X be that portion of the loss 
Q due to radiation. Hence, the ionization loss is Q - X. Considering 
these two loss processes, the probability of energy loss is 



where Wj and are the energy loss distributions for ionization and 
radiation respectively. 

The Landau equation [2], as modified by Blunck and Leisegang, is 
used for the energy loss distribution due to ionization. The distri- 
bution is expressed as a function of Landau's dimensionless parameter 
X as follows: 



The terms used in equations (2) and (3) are defined as follows: 

C^, and x^ are constants given by Blunck and Leisegang [5] and 
are used to fit Landau's distribution to a sum of gaussian functions. 

R is the absorber thickness in cm. 



Q 

W(Q)dQ = f Wj(Q-X)W3(X)dXdQ 



( 1 ) 



0 



where x = + i n -fr - i . 1 1 o 




(3) 
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The quantity "a" is a function of the atomic number Z, the atomic 
weight A, and the density p of the absorber; and b(=v/c) of the 
electrons: 

a = P ' MeV/cm. (4) 

A 6^ 

2 

The quantity b is a correction to Landau's theory given by Blunck 
and Leisegang [5]: 



3.0 

aR 



>A 



In 



m 



2 Ei 



V’ ■ “ 



(5) 



where the summation is over the ionization potentials of the atomic 



electrons, and is the number of electrons with ionization potential 



Q is the average energy loss due to ionization (no radiation) for 
electrons of incident energy E^. , and is given by Sternheimer [11,12, 
and 13] as follows: 



- 



B + 0.43 + In E. - S 



C - a^(X.| - log^Qp/mc) 



m. 



MeV (6) 



2 

where t is the thickness in g/cm , and the constants A , B, C, X, , a , 

^ I W 

and m^ are parameters of the absorber material. These parameters for 
tin, lead, and various other material are listed in reference [13]. The 
parameters for gadolinium are not listed but were determined by the 
following method. 

A^ for gadolinium was obtained by extrapolating from a plot of A^ 
vs. Z/A for the absorber materials listed in reference 13. The result 
was A^ = 0.0624. B for gadolinium was found in a similar manner from a 
B vs. I (ionization potential) plot where I was determined from a Z vs. 

I graph. The result was B = 13.3. C was determined from a semi -log plot 
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2 

of I A/pZ vs. -C, and the result was C = -6.80. The values for a^, 

Xp and nig were more difficult to determine, but fortunately they 

contribute little in determining These constants were numerically 

m 

determined. The values for the term ag(X.|-log-jQ(p/mc)) were 
calculated for tin (Z = 50) and tungsten (Z = 74) in the p/mc range 
(p is momentum, and m is the rest mass of the electron) of the energies 
used in the experiment. Using Z as a basis of interpolation, the cor- 
responding values for gadolinium (Z = 64) were estimated for two 
representative p/mc values (100 and 150). The value for X^ was 
arbitrarily set to be 3.0, and the a and m values were determined 
from a simultaneous algebraic solution based on the interpolated values. 
The results were a^ = 0.418 and m^ = 2.10. The error in the determina- 
tion of ^ using these values was estimated to be }% or less. 

For Wj, the energy loss distribution due to radiation, Blunck and 
Westphal [1] give; 



W3(Q)dQ = BaR(Q/E.)“'^ ^ 



(7) 



where 



a = 1 .4 X 1 0 



■3 pZ 



4/3 in ^ +1/9 



1 



cm 



-1 



( 8 ) 



and B is a normalizing factor = • 

The distribution of total energy loss according to Blunck and 
Westphal is obtained by putting equations (7) and (2) into equation 
(1) and performing the required integration. The result is the energy 
loss distribution for a single electron of incident energy E^. . For 
comparison of theoretical and experimental values, this distribution 
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function was used, with corrections to account for the finite energy 
width of the incident electron beam. This treatment is descriln d in 
Section IV and Appendix C. 
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III. EXPERIMENTAL PROCEDURE 



The LINAC of the Naval Postgraduate School was used to obtain 
electron beam energies from 52 MeV to 92 MeV. This beam was elastically 
scattered at 90° from a thin (3.5 mil) aluminum scattering foil, passed 
through the absorber, and finally analyzed by a 120° magnetic spectro- 
meter described by Kenaston, Luke, and Sones [14]. 

This general experimental arrangement was similar to that used by 
Miller [7,8] and DeLeuil and Raynis [10] with the exception of the 
removal of a 3.5 mil aluminum window and the installation of a ten 
channel coincidence counting system. These changes are discussed below. 

The aluminum window removed was located just before the scattering 
chamber. Since the presence of the window caused a broadening of the 
incident electron beam distribution, its removal decreased the energy 
width of the incident beam for this experiment, as compared to previous 
works done at the NPS. 

The ten channel counting system consists of ten front counters and 
a single backing counter operated in coincidence with the front counters. 
The entire ten channel system has an energy spread of about 3 %. 

The absorbers were positioned approximately 3 cm. from the scattering 

foil as recommended by DeLeuil and Raynis [10]. Absorber thicknesses for 

2 

tin were 1.485, 2.970, 4.455, and 5.940 g/cm ; thicknesses for gadolinium 

2 

were 0.814, 1.610, 3.221, 4.026, and 4.831 g/cm ; absorber thicknesses 
for lead were 2.825 and 4.236 g/cm^. 

A thick aluminum absorber (5.574 g/cm ) was used to correlate the 
results of this experiment to the results of DeLeuil and Raynis [10], 
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which were obtained before the ten channel countin^^ iystem wav . ■ called. 
The results agreed within experimental accuracy. 

Since the electron beam scattered into the absorbing mac.. was 
not monoenergetic, the energy distribution was determined both before 
and after passing through the absorber. Various thicknesses of absorber 
material, including no absorber, were positioned on a remotely controlled 
ladder device. For measurement of the distribution before passing through 
the absorbers (called a zero peak), the ladder was positioned such that 
the beam passed only through the scattering foil. After measuring the 
zero peak spectrum, the various absorbers were moved into the beam to 
determine the energy distribution of the beam after passing through each 
absorber. 

In the case of tin and lead all absorber thicknesses could be measured 
without turning off the accelerator and thus possibly altering the char- 
acter of the beam. This was not possible for gadolinium. Because of its 
cost, only a small quantity was purchased. Hence, only two gadolinium 
absorbers could be run without turning off the accelerator and rearrang- 
ing the gadolinium on the ladder. Since the beam character could change, 
a zero peak measurement was made whenever absorbers were replaced. 

The data represented the number of electrons detected by the coin- 
cidence counting system at the exit of the magnetic spectrometer. A down- 
stream Secondary Emission Monitor was used as a standard for normalization 
purposes, in that each data point corresponded to a given integration 
current, i.e., a certain number of electrons passing through the absorb- 
ing material . 
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IV. TREATMENT OF DATA 



The data reduction for the ten channel counting system is standard 
for multi-channel systems and is on file in the NPS LINAC computer 
library. In this reduction process, three characteristics of the 
counting system are required. The characteristics are: 

(1) the energy spread of the counting system, 

(2) the energy seen by each front counter*, 

(3) the relative efficiencies of the front counters. 

With this information the actual energy corresponding to the front counter 
data is calculated. Counting rate and background corrections are also 
performed during this computerized process. 

In the lower energy regions of the spectrum, where all counter read- 
ings should be nearly equal, a counting discrepancy was noted in that the 
lower energy front counters read consistently higher than those on the 
higher energy portion of the counting spectrum. This effect occurred at 
energies lower than required to obtain distribution half-widths and thus 
did not affect the data reported. However, it made confirmation of the 
shape of the distributions at very low energy {< 20 MeV) impossible. 

The electron beam energy is not monoenergetic as required by Blunck 
and Westphal theory. To compare experimental and theoretical results, 
the zero peak data must be used to modify the theoretical predictions. 

With no data treatment, the theory will predict half-widths which are 
too small and, in the case of asymmetric zero peaks, erroneous most 
probable energy losses. Removal of the aluminum window noted in Section 
III narrowed the zero peak to the extent that it could be treated as a 
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symmetrical gaussian distribution. The IBM 360 computer of the Ni'S 
was used to unfold the zero peak energy distribution into the iheoV'y. 

A histogram method described in Appendix C was used for this unrolding. 

The computer program used to accomplish the unfolding is Appendix D. 

In the comparison of experiment with theory, the measured incident 
energy distribution was unfolded as previously described, properly 
normalized to the experimental data, and plotted. From these plots 
the theoretical half-widths and most probable energy losses were 
determined. The experimental data were plotted along with the correspond- 
ing theoretical curve, and measurable parameters were compared. The curves 
are shown in figures 1 through 38 of Appendix B. 
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V. RESULTS AND OBSERVATIONS - 



The values predicted by theory and the experimental results for the 

most probable energy loss and half-width for tin, gadolinium, and lead 

are shown in Tables I, II, and III respectively. In these tables, BW 

refers to Blunck and Westphal theoretical predictions, Qp represents 

the most probable energy loss, and HW is the half-width. Data are given 

for predictions with and without beam folding. 

On comparing the theoretical predictions with beam folding and the 

experimental results, it is seen that good agreement exists for target 

2 

thicknesses less than about 3.0 g/cm for both half-width and most 

probable energy loss. At thicknesses greater than this, the agreement 

in most probable energy loss is still reasonably good, but there is poor 

agreement in half-width. The greatest variation in most probable energy 

loss is 8 % while the average variation is 3 %. These percentages, while 

averaged over all material, are typical and not material dependent. The 

average half-width variation is 9 % where half-widths were obtained. As 

can be seen from Appendix B, figures 2, 4, and 6, the theory does not 

predict a half-width for a certain target thickness for each element. 

This cutoff thickness is smaller as the atomic number increases. For 

2 

example, the theory would predict a half-width for 5.6 g/cm aluminum, 

2 

but not for 4.2 g/cm lead. 

It is concluded from these results that the Blunck and Westphal theory 
predicts correctly the most probable energy loss for all targets used and 
predicts satisfactorily the half-widths for the thinner targets. Thin 
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targets can be defined empirically as those satisfying the relation 

*] yo 

T(Z) ' $ 13.5, where T is the absorber thickness (g/cm^) and 1 is 

the absorber atomic number. 
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TABLE I. Energy Loss Distribution Characteristics of Ti 
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TABLE II. Energy Loss Distribution Characteristics of Gadolinium 
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** Effective thickness used due to 8° error in target angle. 



TABLE III. Energy Loss Distribution Characteristics of Lead 
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APPENDIX B. Figures 
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Figure 
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Figure 





CTj 

u:; 

LO 



Q 

CTl 

cj:j 

XT 



ezi 



o 

uT^, 

rn 



CIj 

o 

CP 

rM 



cn 

o 



Q 

a 



e_j 

p- 



Cj 

di 



L4.J 



Sj.MHCD □ J^ribMldOM 



35 



Figure 10 
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Figure 11 
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Figure 12 
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Figure 15 
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Figure 18 
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Figure 20 
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Figure 21 
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Figure 25 
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Figure 27 
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Figure 30 
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Figure 31 



C3 

CJ 




CM 

—I 



a 

a 

CO 

cn 



a 

a 

TT 

CD 



CD 

a 

r 

CD 

rv 



CD 

CD 



CD 

in 



3 

UJ 



CD 

CD 

CXI 



CD 

CD 

OO 

CM 



CD 

CD 

XI* 



CD 

CD 



Siwnoj oizndkiyoM 



57 



Figure 32 
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Figure 33 
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Figure 34 
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Figure 35 
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Figure 37 
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Figure 38 



APPENDIX C - BEAM FOLDING TECHNIQUE 



The theory of Blunck and Westphal [1] assumes that the beam of 
electrons striking the absorber is monoenergetic. The beam of electrons 
produced by the NPS LINAC, or any other linear accelerator, for that 
matter, has a finite energy spread about a maximum or most probable 
energy point. The fact that monoenergetic electrons are not avail- 
able to strike the absorber must be taken into account in computing 
theoretical predictions if a meaningful comparison with experimental 
results is to be made. This has been done here by a technique termed 
"beam folding". 

Beam folding is accomplished by a number of well defined steps. 

Energy distribution curves may be approximated by histograms. These 
histograms consist of a series of "bins" of area W(Q)aQ, where Q is 
the energy loss and W(Q)aQ is the probability of loss between Q and 
Q + aQ. Thus, each bin has an "address", Q, on an energy coordinate 
scale. To accomplish beam folding, the beam distribution must be known. 
This is observed experimentally and approximated in the computer by a 
gaussian curve of appropriate half-width. The energy at which the 
maximum occurs is established by the energy of the beam and the magnitude 
of energy loss incurred by elastic scattering of the beam as it impinges 
on the thin aluminum scattering foil prior to striking the absorber. 

The width AQ is then selected for the predicted distribution. This must 
be a small, but finite number where numerical techniques are to be used. 
This same width is used to break the beam distribution into a histogram. 
This is illustrated in figure 1(a) of this appendix. The reason the 
same width is used for both distributions is to facilitate computer 
programming. 
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(c) Energy Distribution of Eg 
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Each bin of the beam distribution is now treated as a monoenergetic 
beam with energy commensurate with its center location on the energy 
scale. Each bin also has a definite magnitude, or weight, with the 
magnitude of the center bin being unity. The formulae of Blunck and 
Westphal is now applied to each of these beams and an absorber distri- 
bution curve results for each, with a maximum amplitude proportional 
to the height of the appropriate beam distribution histogram. Each of 
the absorber distribution curves thus obtained may be thought of as 
being plotted and added to previously determined curves, using the 
energy of the electrons as a correlation point, as depicted in figures 
1 (b), (c), and (d) of this appendix. This is accomplished on the 
computer by adding the contents of each bin of the same address and 
plotting the cumulative total. The bin of address E is depicted in 

A 

figures 1 (b), (c), and (d). Note that the beam distribution for histo- 
gram E4 is centered over the maximum point of the distribution. This 
is important as a false picture could easily be presented if the histo- 
gram were not symmetrical as the distribution would then appear skewed. 

Since the half-width for a particular absorber increases as the 
thickness (g/cm ) increases, the beam folding technique is affected in 
that thick target energy distribution, whose half-width is large, will 
be affected much less by the finite beam distribution than will the 
energy distribution of a thin target whose predicted half -width is not 
more than 2 or 3 times the size of the beam distribution half-width. 
Thus, in selecting bin width (aQ), it is desirable to make the width 
proportional to the target thickness so as to apply the beam folding 
technique in greater detail to thin as opposed to thick targets. The 
computer program (Appendix D) applies the target thickness in direct 
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proportion to establish a bin width for subsequent beam folding 
calculations. Thus aQ is much smaller for the thinner absorbers 
compared to thick ones. 
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APPENDIX D. Computer Program 



C PU'JPnS':: THT*^ OPOGPAM ACCFPTG PAc-AVPTFpG FPOM ckjtpqy 1 DSG 
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C WILL PLOT BOTH THF F XP PP I »BEMT AL ANO THTorETICAL cm^pGY 
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r. 

C SCnpF; THE PROGRAM IS SET IIP Tn MAKE F, OMPUT A T I pm R qc|_/vtiv = 
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C AND LEAO^ 

C 

C METHOn: THP PROGRAM PL'^TR A Theory FIIRVP BASEO org TH^ ^1^ 

C OF GAUSSIAN GURV^S PPRFPRMcr) ry BLUNCK AMO WESTPHAL '’PR 
C THE PREOTCTPO ENjcRny niSTOIOUTIlN CURVES. THE INTcgRA,- 
C TION OVER EN^doy LOSS INOICATEO TM THF A Epp F Mc mt i OM<= 0 

C PAPf^R IS PERFORMPP ON the COMPUTER OSTMG POINT SAUSS 
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C B-W PR'^DICTIO^I WITH A MO NE N= R 0,1= T f C, pc PL'^CTRONS A\ip 

C THE OTHER, WHICH ALSO ‘JAS EXPERIMENTAL RESULTS thcrcom^ 

C represents the same TKcpRY WITH A PO LYE me R GE T I C n<^AM 
C DISTRIBUTION FOLDED IN., 

G 

C THIS PROGRAM IS WRITTEN USING FORTRAN IV I ANGIJAGE AMD tmc 
C NPS COMPUTER FACILITY PLOTTING PACKAGE, TOTAI RUNNING 
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lWUF{BOP),XUF(BPP),EXX(RnO),EXY( 2 PP),IXYI 2 CO),TR( 4 ) 

EXTERNAL FCT, BEAM, WI ON , WRAD 



D01200 J=ltF^ 
E( J )=0 
WF{ J)=0 

001201 J=l»2O0 
EXXI J) =0 
EXY{ J)=0 
I XY( J) =0 
001202 j=i » ^on 

WOT{ J) =0 



HW( J ) = P 
XXI J )=0 



WUFI J) =0 
120? XUF(J)=0 



READ IN parameters FOR EXPERIMENT, 

FR beam energy 

7 ATOMir number 

A ATOMIC weight 

T target THICKNESS: GM/SO C.m 

07EP0 SPECTPOMFTFR INDICATION PE LOCATION np rcam 

DISTRIBUTION ( Mcy) 

CTSNQR number Tp WHICH 'EXPERIMENTAL COUNTS WILL BP 
normal I 7 En 

HT HALFWIDTU OF REAM DISTPIRUTIPN 
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<?'=A0(5,in FR, A,T,0ZE'?0,CTSMnR,HT 

HTVIOE 7 <=Rn PFAK INTO ••RIMS" IP WI OTH O^LX, WIRTH 

INCREASES with tapGFT thICKNFSS, 

n‘"LX =.0 3’*'T 

IF( T,LTo 1. ) nELX = . O^'-T 

OAVF IS THF RFCniL LOSS np Pt'^CTRIMS STRIXIMG AIU*'’IM')m 
NUCIEII OF THP SCATTPRTMG f^IL. I0NI7ATION inss IV 
SCATTERPP IS NEGLECTFO, 

0AVF=FR^F0/ (93i , 4 <-a » 

C'^PR IS THF rnRRFCTIOM TT Rc A9n‘=n TO' EXPERIVpMTAL FMCi^fSY 
TO rORRFCT FOR S P'' CT ROMF TF R frpor^ 

CnRR=07.FRn-QAVF 

12 F0RMAT(5(lX,Fo,3,|y,x6) ) 

PFAO in up to 2'^n FXDfrimfmTAI ROIMTS; fxX is FMFDOy 

AND IXY IS CntJMTS, THERE MOST RE 4.3 OATA CAROS AVO tmp 
LAST EXX MUST RE ZF^n, 

RE AO (5, 12) (fxx( J) ,IXY(J),J=1,P00) 

ESTARLISH PARAMETERS FOP SCATTERING FOIL-, 

FOIL = ,0'»4 

7FOIL=l R, 

KSIG=0 

no I J=1,20P 
1 EXY(J)=IXY(J) 

M = 0 



El IS value of beam energy INCIOENT on target, 

FI=fb-OAVE 

TEMP=EI 



THE following PORTION OF THF PROGRAM UNFOLOR REAM otstpT- 

bution, icent is am integer valuf assigned to tmp 

HISTOGRAM CENTfrED OVFp THE PFAX OF THE DISTpIp.U- 

TIOM, EACH HISTOGRAM IS DFNOTPD py ITS fnERGY p AND 
MAGNITUDE RELATIVE TO TH^ rpaK OF THE DISTRIBUTION, WF, 

ICFNT= ( HT+OELX )/OELX 

X = EI-TCFNT^nF|_x 

202 X=X+DELX 

IF(X,GT„ (EI+HT) ) GOTO 1003 
M = M+1 
E(M)=X 

FUNCTION BFAM RFPLICATfS THF D T S Tp i bhT I ON or RFAm F\|cpf;v 
AS A SIMOLF GAUSSIAN CURVE BASFO ON THF Ex pp R T mf NT AL ( Y 
OBSFPVFD HALFWIDTM AND IS LOCAT'^D PFLATIVF TO THF Ri^AM 
ENERGY BY THE COMPUTFO AVFRAGE RECOIL LOSS, 

WF(M) = BFAM( EI,HT,X ) 

GOT0202 

1003 WRITF(X),?2) OAVE,HT,fou ,7Ff.lIL ,fi 

no loof L=i ,M 
nx=o 
N = 0 

lOOF N=N+1 

c 

C USING faCH histogram AS A "Bfam'^, THF FOLLOWING POPTIOm of 
D THF program COMPUTFS an fncrgy niSTRIBUTIOM pop EACH 
C OF these "RFA*^S'^ and STORfs THFM IN "BINS” LARFLED py 
C ENFPGY, XX, THF rUMlILATIVF AMOUNT OF rlfCtrom COUNTS IN 

C THE BIN IS WOT, 0 IS thF FNFRGY LOSS OVER WHICH TMF B-W 

C PPEOICTION IS INTFGRATED AMD DX IS THF amount BY WHICH 
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0 CMANGF*; FDR EACH SI)FCF<;s TV T ^ 1 T AT I OM, IMT'^G°ATinM 
IS CAOPIFf) OUT IN SIIBROUTIMP 0 QG 3 '’> 

DX=DX+nFLX 

0 = 3 ., O^T”^ ( 7 / I 3 , ) V^-OX 

IP( 0 »L‘^oEI IGOTOIOOA 
N = 0 

GOT01005 

lOOA. IF(0,LT,o)r,nTnior>6 
EI=P(L) 

XL = , CO'^000'^1 ^0 
VALI M=,V‘^'«^0 

CALL n 0 G 3 ?( 0 -,,XL , <=C T, 7 A , F I t 0 , ^ , A , Z ) 

CALL nOG 33 ( XL, VALI M, Ff T, 7 B, P I , 0 , T, A , 7 1 
CALL OOG^?(VAl.rSO , PC T , ZC , E I , 0 , T , A , 7 ) 

WO= 7 A+Z 8 + 7 C 

W 0 T( L + N-T ) =W 0 T n +N-I ) + W 0 *W F ( L ) 

XXTL+N-l )=F(L)-0 
J=L+N-1 

the 0TSTRI8UTI0N pno jhp HISTOGPA^ CORR pR PONO T NG Tn THP 
CENTPR OP thf 8FA*'^ 0 1 STP I8UT I0^' IS ppmovEO ANO STORpo In 
WUF ANH XIIF AS THE UNPOLHEO COUNTS AND PNPRGY 
RESPECTIVELY, 



I F( L ,NP., ICFNT )GQT01005 
WUF( M) =W0 
XUF<N»=PI-0 
K = N 

GOTO isns 

1006 CONTINUE 

on !717 1=1,7 

ELIMINATE LOW NUMRpRPO "PINS" AS THPY ARE ONLY PARTLY FTI- 
LEO ANO ARE NOT REPP FSpNTATI VP OF THE CU^ULATIVp 
DISTRIBUTION, 

WUF( I)=WUF(8) 

171'^ WOTI I > =W0T (8) 
on 6 1=1, A 
6 T3( I )=T6( I ) 

GOTO 77Q 

777 00778 N=1,K 

PROCESS THP UNPOLOPn tHPORY CURVE, SPTTTMg EXXCM pQUAL 
TO ZERO ALLOWS PLOTTING OF THIS THEORY CURVE WITHOUT ANY 
EXPERIMENTAL DATA, 

WOT{N)=WUF(N) 

XX(N)=XUF( N) 

KSIG=1 
J = K 

778 CONTINUE 
EXX( 1) =0 

no 5 1 = 1, A 

5 T 0 ( I )=T 6 ( I ) 

779 WQTMX=0 
001C07 M=3,J 

OFTPR'-’INP THP MAXt^'UM POINT OP ^hc 01 S TR I RUT ION, WQTMX, 

AND THUS PSTARLIS'-' THp *^0ST PROBABLE fN'pRGY LOSS, QP, 

IF ( WOT ( M ), I, P, WQ 7 (N -1 ) o ANO, WOT( N-i ) , GE . WOT ( N- ? ) » 
IGOTOIOOB 
GOTO lOO"^ 

1003 IF (W 0 T(N -1 ) ,|.T, WOTYXIGOTO lO-'T 
WOTMX=WOT ( M- 1 ) 

FI =TFMp 
OP=P T-XX(N-i ) 

1007 CONTINUE 

NORMALIZE ALL PniNTS OF OISTPIBUTTHN ANO STHRF TM HW, 
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DO inoc M= 1 , j 
HW{M)=WOT( M)/WQTVX 
CONTIMUF 
J = J-T 
HWI=0 

C DFTFRMJMC HALFWrnTH np jhc TH^HRcyifAL HI STR fRUT I HN), HWA . 
DO imo M=1,J 

IF{ HV-IC GT. •, t;-, AMD, HW( M), |T, ,5)HWl=XX(M)+( ,F- 

IHW(M) ) « nFLX/(HW( ■^+1 )-HW( V) ) 

IF(M,PO, J„ AND,HWI, '^0>n JOr'TOl iq 
GOTOl?! 

11® WRITF{6,?0) 

HWA = '^ 

GOTO 101^ 

20 format (/,? X ,• NO HAL'"WIDTH OBTAIMfOM 
1 ?1 I (HW{ M ) oGT,>, 5 > AMD,, HW( M + 1 ) , L T > , F ) HW A = X V ( M ) + ( H VM 'M 
l-„ 5) ’‘DTL X/ ( HW( M) -H W( ^^+1) )-HWT 

1010 rONTTMU^ 

ESTABLISH CUTOFF SIGNAL ^OR PLOTTING SUBPOUTINP. 

XX( J+1 1 

ESTABLISH CORRFCT GRAPH TTTLFS FOR ABS<^<^BrR USFD. 

IF(7<,GT,13,1)G0T02R 
DO 0 1=1 ,4 
3 TR( I )=T3{ I ) 

GOTO 1011 

?P lF(ZoGT,2P,l 1GOTO50 
DO 101=1 ,4 

10 T?( I ) = T1 0( I ) 

GOTOlOll 

50 1F(7„GT,50,1 )G0Tn64 
DO 8 1=1,4 

8 T2( I ) = T8 ( I 1 
GOTOlOll 

64 TP( Z^GT, AA, 1 )GOT082 
DO 9 1 =1 ,4 

9 T2 ( I )=T9 ( I ) 

GOTOlOll 

82 DO 7 1=1,4 
7 T2(I >=T7 ( I ) 

send DATA TO PLOTTING SUBROUTINE FOR PROCESSING, 

1011 CALL GR APH( XX, WQT, FXX, <=XY, T4, T3, T2,T1 ,HWA , FI , T, 
1C0PR,CTSN0R,QP ) 

PRINT RESULTS OF COMPUTATIONS IN TARULA^ FORM, 

WPTTF(6,2?) Qd,HWA,T, 7,TI 
IF(KSTG,f^0,O) GOTO -^77 

11 FORMAT{7f^l0„0) 

27 FORMAT!/, 2X, ' QO : ',F7^4,' HAL^ WIDTH: ',F6,^, 

1* T(GM/SO-C‘M : «,F6-,4,« 7: •,F^,0,‘ INCIDENT EM'^RGY: ‘ 

2,06. 2,/) 

STOP 
END 
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suBPOijT INP nor-3?(xL,xij,FCT, Y,c } ^n,T,n,7) 



PIJPPOSF: SUBROIjtimf pppe^pMS THP IMIpb^MIOKJ p-QiiIPph 

IN THF R-W THf^TRY, 

IMPLICIT RFAL^-RI ) 

A = , POO" ( XU + XL ) 

B=XU-XL 

O-^CPP^TAOTRnO* R 

1 ( FCT ((A+C),7,n,T,FI,0)+FCrn&-D,Z,n,T,'^I,0)) 
C=,4P?«n575S77pf^-c.i 700^-B 

Y = Y+,B1771P‘'7fB4F?'^3Fn-7’' 

1 (FCT((A4.C) ,7,D,T,PI,6)+PCTnA-C» ,Z,n,T,pi,o) ) 

C = ,4B2’B!17 77q77 53?7n0-^B 

Y = Y-*-o 1 ?6PiSn77f,5/i A,7 ] o^oO-l 

1 ( PCT ((A+C)t7,n,T,ci,0)+PCT((A-r),7,0,T,FT,0)) 

C = , 4674 5 3n7 796 0R6°fl4nn «R 
Y=Y + „ 171 3 BO?i4PA'^ 107170-1* 
l(FCT((A+n,7,0,T,PT,0»*FCT((A-C>,Z,O,T,^I,Q)) 
r. = „ A4R1 6OS77RR7O76O4O0"'8 

Y = Y+,?1417O4oniln 334r0-’'- 

l(FCT((A + C»,7,n,T,FI,0)+FCT((A-C),Z»n,T,c:i,0)) 

C=<,4 24^3 3B''f R6'S7R4qQ0n*B 

Y = Y+o2'^4oqq7q631 ’ RRoaRO-l '* 
](FCT((A+C),7,D,T,PI,Q)+FCT((A-C)»Z,n,T,t:i,0)) 

C = „ 397741 BQ79R7071 ponofeR 

Y = Y+ , 79 3 42 0 467 3^2^7 7 740-1 *r 

1 (FCT((A+C) ,7,n,T,Fl,0)+FCT((A-C),7,0,T,cy^g)) 
C=«3660P105P37Ol 44B400-«:R 
Y=Y+a37qiini3B«ie0 9230-l* 
l(FCT((A + C),7,0,T,PI,Ql+FCTUA-C),Z,0,T,FI,0n 
C = -» 33 152 21 334 651 0760O0*R 

Y = Y + ^36177P97 0F4474253 0-1»- 

l(FCT((A+C).7,n,T,FI,Q)+FCT((A-C),Z,O,T,PI,0)1 
c=, 203 R57B7»62 0351 1 600 *R 

Y = Y + o 390 06 047 807 53 5 157 9-1 4- 

I (FCT((A+C) ♦7,D,T,FI,0)*FCT((A-C),7,0,t,ct,q)) 

C = , 2 534499 5446 61 14 7C00-^R 

Y = Y + »41 655 0621 1 34;73 7780-1 * 

1 (FCT ((A + C),Z,P,T,Fi,Q)+PCT((a-C),Ztn,T,'^I,0)) 

C=, 2 106 75 67 «0 657 176700 -6p 
Y=Y + «43 876 04A 607 20 1906 0-1* 

1 (FCT ((A+C),Z,n,T,Fl,0)+FCT((A-C),7,D,T,EI,Q)) 

C = . 1 65074701141 '7638700*8 
Y=Y + « 45 5 860 393 47 88 1047 0-1 * 
KFCTMAi-CK^.n^TtOi^Qi+ccTlCA-OtZtOtTtOItO)! 

C=, 1 1964368 11 260685400*8 

Y = Y + »46077 1 90 5^04 0 7 28 3 0-1 * 

l(FCT((A+r) ,7,0,-r,EItO)6FCT((A-C) ,Z,0,T,FI,0) ) 

C=*7 7235 98 079139825 0-1 *8 

Y = Y+, 47 Ri 9 3 60 07 0 63 7 430 0-1 * 

1 (FCT ( ( A+C» , Z ,0,T,F I ,OH-FCT ( ( A-C ),Z»0,T,fi,q)) 

C = , 2 41 5 3837 P43 869] 570-'' *P 
Y=R*(Y+„48?700447577679ao0-1* 

1 (FCT ( ( A+C) , Z ,0,T, E I ,0 ) +FCT ( ( A-C ) ,Z ,0,T,F I , Q ) )) 

R'^TURN 

END 
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RPAL FONCTION W I ON ^ ( X , I , T , 7 , A , 0 ) 

PURPOSE: THIS •^UMCTIOM CALCULATFS THAT PART nr thp PNFpr.Y 
0 ISTR I RUT T ON OF AH P| ‘^rTRPU OIJF To L^ss nr fmcp^y «y 
IOMIZATIOH WHTLF PASSTHG through an ARSORRcp , 

IMPLICIT RF AL<^R{ A-H,0-Z) 

TL ARE the average lOMIZATinri pnTFMTIALS P<=R eleCTROM FOR 
LEAD FOR FACH SHELL, BEGINNING WITH |< rhclL-. 

NPOTL IS THE MU'^BFR OF IONIZATION POTENTIALS USED TOP |paD 

IN THIS CONTEXT, L = LEAD, G = GADOl I N J !JM, S = TTN, C = COe>opp, 

AND A=ALUMINU*S 

REAL*p IL( 6 )/^ OBRono,„ 01 470n , , , , onr.490, 

U000076 , . 0^0002/ 

REAL ’'ft NLI G)/2, , 9-. , 1R„ , 1 B-, , / 

REAL-'^B IG(5)/,0SOR4,,nn7B5,n001o9,,0DC?4,,00 007^/ 
real * 8 NGI S) /2, , Bo , \B, ,o, / 

REAL'^=8 T SI E )/, 07 02, ^ on A4 , . OOP 7 7, o 00 006 2, , 0 conn 1 / 

REALtR NSI S ) /2o , ,19o ,1 R . , Ao/ 

REAL-^B IC( ^ )/,0OB9B,, ooooo, ,OOOn6B/ 

PEAL*8 NC ( 8 )/2, , Bo , IR, / 

REAL*B I A( 2 ) /, 00] 86, , 00Cn8«/ 
real '< 3 NA( 2 ) /2o , B, / 

NP0TL=6 

npotg=r 

NPOTS=e 

NP0TC=8 

NP0TA=2 

SB IS THE SUMMATION OF IONIZATION PHTEHtiaLS eQR THE 
FXPERIneNTAL ATDMo 

SR=0, 

XMASE is the rest mass OE an electron in npv, 

XMASe = o 511006 

BETA IS THE NORMAL V/C USEO IN RELATIVISTIC CALCULATIONS, 

GAMMA IS THE NORMAL TERM USED IN RELATIVISTIC CALCULATIONS 

8ETA=0S0RT ( ( ET ^-P I + 2o *:E P XMASE) / I E I -r f I +2 , ^ f jo x MA S F 
1+XMASF^XMASE) ) 

GAMMA=1,/( lo-BETA* «-?)**( 1, /2, ) 

AT THIS POINT, THE PROCRAM MUST CONpmTf THF proper rON- 
STANTS FOR THE FXP FP I MFN’^AL 7, A SERIFS OE LOGIC STATF- 
MENTS select TH^ CORRECT FORMULAE, 

13 IF( 7.GT. 13, 1 IG0T02O 

compute SB FOR FXPrp Jmf^jtaL Z, R IS THE IONIZATION POTEN- 
TIAL from ome SHELI.» 

00 71 T=l,NPnTA 

R=IAm^NAn)*nLnG('>,-^EI/(IA(n’^-(l,-8ETA’*-'^2))) 

71 SB=SB+B 

ESTABLISH VALUES FOR R2 AND A], THESE ARE MATfpjai OE- 
PENOFNT constants USeo TO CnMOUTr THE AVERAGE FMeroy 
LOSS BY TONIZATIDN ANO ARE SUPPLIEp BY the RTEPNMriMFr, 
PARER, these ape R ANO A, RESPECTIVELY, FROM HIE pAPER, 

R2=J 6,T7 

A1 =.07-i0 

PELT IS A PART OF THE EXPRESSION BY FTEONHEyvcp FOR THE 
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AVFPAGE FNJCRoy l.nss Dlj'= TT inMl/ATTOM, NIJMRFPS HcprjN 
PPPRFSFNT OTHPR s T ER 1 mf o CONSTANTS AS FOLLOWS: 

small a „t> 00 (S 

SMALL M 3,F10 
XI 3.000 

-C 4,P10 

DELT= 4 , 4 O 6 -<^nL 0 Gi 0 ( RCTA + GAMMA)- 4,21 f-, 000 G^ 

1 ( 3,-DLOGlO ( RFTA *- gamma ) ) ^^3 FI 

GO TO in 

2 ° IF(7 »GT.no , IIGOTOFO 
no 72 1 = 1 , NPOTC 

fl = ic m*MC n )^DLnG(R. IC( I ) ' (1,-R‘-TA^*^'») ) ) 

72 SR=Sn+P 
B?=l 5 ,oq 
Al = .07'^l 

DELT =4, f,06'^0| OGl 0 ( BETA -^G AMMA ) - 4 , 74+,, n PO* 

1 (3o-0L0G1.0(PFTA-^GAMMA) »-^‘'3>7«. 

GO TO in 

50 TF(7oGT,F 0, nG0TO64 
no 73 1 = 1 , NPOTS 

P=IS ( T »*MS ( I )^DL0G( p,^fi/( IS(n'^(l,-RFTA^'^?n» 

73 S3=SR+B 
B3 = l 3„ 93 

Al =, 0647 

OELT =4 -. 606 ''OLOGIO( BETA^GAMMA 1 - 6 , 29 + 3404 <'( 7 ,- 
10 L 0 G! 0 ( BETA’ GAmma) ) 

52 

GO TO in 

64 IF(7 ,GT,64, nGQTOPP 
on 74 1=1 , NPnTG 

B= TG(I » "'MG ( I )*DL0G (’nFi/ciGdl^'M ,-BFTa*'^2 » ) ) 

74 SB = SP-*-B 
B2=13,3 
Al=a0674 

0 FLT= 4 -, 606 ''OLOGI 0 ( BFTA^GAMMA )- 6 , B +,41 «*( p . - 
IDLOGIOC BETA'- GAMma) ) 

n^' 2 * 10 

GO TO 111 

B 2 TF( 7 .GT, 92 .,nGnTnBAo, 
no 75 1 = 1 , MPOTL 

R=IL(T )TNL( I )* 0 L 0 G( 7 .,'^EI/( IL (1 )-^(] ^-BETA^-' ^l ) ) 

75 SB=SB+B 
B 2 =l ?• 91 
A 1 =,0609 

0 FLT= 4 o 6 06 *nL 0 Gl 0 ( RET A ^^g AMM A ) -6 . Q 3 0 65 2 -f- 

1 ( 4 »-DL 0 G 10 ( PFTA’f GAMMA) )^’^P, 41 
GO TO 111 



AR IS THP OHAMTTTY SMALI A, FROM PIONCK AMO W^STPHAL, 
MULTIPLIED BY R, THE TARGET THICKNFSS IM C E N T I ME T^" R S , 

111 AR=0, 154*7 *'T/A/BETA/BFTA 

COMPUTE BS7, the small P SOUAREO FROM THE BLUNC!< AMD 

lfis^gamg paper. 



BS2=3.^SB7 ( 7*AP) 

VARl , VAR7, AMD VARB AR^ COMVFMIFNCF STORAGE LOCATIOMS 

USFD TO STORE COMPUTED PORTIONS OF 5 T <^R MHF I M ER ' S AV^RAC^ 
ENERGY LOSS, 

VARl=AnT/RFTA -*"7 

VAR7=B24,4P + R, fcPLOGI RETA-^GAMMA ) 

VAR3=DL0G( ^I )-RFTA**p 

QAVF IS THF FINAL VALUE FOR TH^ AVERAGE ENERGY LOSS. 

0AVF = VAR1* ( VAR2*VAR^-DFL'^ ) 

Bl=DSQRT(BS2) 

C 

C X IS THF PORTION OF SOM^ TOTAL ENERGY LOSS 0 DUO T'^ 
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PEAL FtjNCTION VIR AD^ ( X , p I , T , 7 , A , 0 » 

PURPOSE: THI5; FlINrilOO OALCIILATFS that part nc tmF '-mtogy 
0 ISTR I PUT I niM np aU FIPOTR'^U DI.IF T»^ LOSS 'iF pnFROY PY 
RADIATION) WHILF PASSTNC7 T||OOUOH AN) AHSORRi'P. 

P^PLICIT REAl -NPI A-H,n-7 I 

U IS AN INTFRMEDIATF STORAGE LOCATION USFD IM C^MPUTIMr, 
ALPHA, 

U=lP3,0/Z*-"( 1. /T,0 > 

CnuouTc VALUE OP ALPHA ‘-p FOR TARGET THTCKN)pSS. ATDN'IC 

NUMPPR Aun ATpN'ir WFIOHT of cypno I yFMT aI ABROprcr, p IS 
THE TARGPT THICKNESS EXPRESSED IN CFNTIMPTFRS. 

ALPHR='^,001 a-nT’'-Z ''7/ A''( A, 0/?. , ON olOG(U)+1 

compute VALUE FOR R, A NIOPMALIZING Cp»NSTANT, 

R = l-, /DGAMNA(ALPHP+1, ) 

COMPUTF WRAD, THP mumper HF COUNTS ^XPECTEn AS A PPRUI ^ OF 
RADIATIDN FNcpf,Y I'^SSES, X RPPRES''NTS THAT PORTin^ n>- 
SOMF TOTAL ENERGY LOSS 0 WHICH IS NOT LOST BY IONIZATION 

WRAO = (R*ALPHR'N(x-n-^''ALPHD) / ( Et’^>! alphp n/x 

RETURN 

END 



REAL function ROAM' R( F I , HW, X ) 
implicit R FAL*p ( A-H,0-Z ) 

Y = P, 

ALFA=4, tOLOGI Y I/inW^^HWI 
BEAM=DFXP(- ( IX-FII nalFA) 
RETURN 
END 



R‘=AL function FCT^pIX, 7, a, T,cI,0) 

pilopnSF; THIS FUNCTION! Mi)L T l PL I C AT T VF | Y JOIMR tpGFTH'’R th = 
prediction of fncpoy loss D'JF TO IONIZATION AND THF l orf; 
DUE TQ radiation, 

implicit RFAI *p( A-H,0-Z) 

FCT=WRA0( X ,F I ,T, 7, A , 0 I ■'•wi ON (X , FI , T , z , A, 0 ) 

RFTUPN 

END 
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SURROUT INF GRAPH ( TH X , T HV , F X X , X Y , T 1 , T ^ , T4 , MW , T , T , 

ICORP.CTSNOPfQP) 

C 

C PIJPPOS'^: THIS SURPnUTIN^ PL<^TS THF RFSMLTS Oc thfhdy CA|_- 
C rULATIONS AND F X P I MFNT AL Rf^SULTS. 

C 

IMPLICIT R cal*P( A-H, n-7 > 

RFAl -^P TITLi^XI 1 ) /• MFV • / 

R = AL'*'R TTTL'‘Y(pI /'MORMALIZF n COUMTS '/ 

RFAL-*^B TITLFOI P) /• nnx IP, PGV L'PSS'/ 

OIMFNS THN TI ( A) ,T? ( A) , TR( 4 ) , T4( /V ) , thX( ■>0^) ,THYI p'^0) , 
IFXX(POQ),PXY(?00),CTSUP(?pO),CTSDN(pO^) 



DFTERMINF SCALING FACTORS AMD CHANGE DATA TO OISplACFYFMT 



C 

c 



PO 

91 

20 

21 

22 

2 ?? 

100 

30 

31 

32 



3R 



AO 

101 



FO 



X = 3, 

Y = F« 

STXIN=( F I-QP)T?> 

ISIXlN = SIXIN + o>=: 

SIXIN?= ISI XIN 
SI XI N= SI XI N?/2, 

IF(HW«FO,n)GOTnoo 
XNCR2 = , 5-^HW 
INCRP=XNCR P + ^ 5 
IF( INC.R2«LTo1 ) IMCR? = 1 
XINCR2=INCR2 
XINCR=XINCR?/2 . 

SP = SIXIN-6,'^XINCR 
IF( SP.»LT,0)GnT0P0 
GQTngi 
$P=0 

TNCP = < FI-1, 5*T)/^o + ,5 
XINCR=INCP 
N=0 
1=0 
1 = 1 + 1 

IF(EXX( I ),GT,1, I GnTO?l 
GOTO?? 

N=N+! 

GOTH PP 

IF(N,F0,O)GnT0100 
on 23 J=1,N 

FXX( JI= ( FXX ( J ) -SP) /X INCR+CORR/XINCR 
IF(J*F0,M) GQTQP3 
CONT INHE 
1=0 
M=0 
1 = 1 + 1 

IF(THX( n* GT, 1, ) GHTO 31 

GnT032 

M=M+T 



G0Tn30 
DO 33 J=l,M 



THX( J) =(THX(J)-SP) /XINCR 
IF( J,FO,M)GCTn33 

IF(THY( J+1 ),GT.THY(J) ) THYM A X= T HY (J + 1 ) 
CDNT INUF 

IF(N«EO,n) GOTOIO! 



SCALP=A, 

IF(T»GT,A, P) SCALF=P, 
no 40 J=1,N 

FXY( J) =FXY U )/CTSNGR + SCALF 

CONTINUE 

DO 50 J=I , M 

SCALF=A, 

IF(T,GT,A,P) SCAI,F = ^. 

THY( J) =THY(J) /THYM AX ‘'RFAL E 
CONTINUE 



SORT EXX AND flimimATF VALUES OFF GRAPH 
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\ 



IF(NoEO,^) GOTnlO:’ 

K = N-1 

no so 1 = 1 ,K 

I p 1 = I + 1 

on 6 0 J = IP! ,N 

I''(FXX( T )> I PXX ( J ) ) goto 6"^ 

TE'1P=FXXn ) 

T<Em=FXY ( I ) 

Exx( i) = '^xx( .n 

PXY( T) =f^XY( J) 

FXX( J)=TC^MP 
EXY< J) =T^:m 
^0 CONTIMUF 
J = 0 
KK=0 
K = N 

m 61 1=1, K 
IF(KK^F0„1 IGHTP^l 
I f C =XX ( I »o LToO )nnT061 
J=J+1 
N = J 

EXX( J)=FXX( T ) 

= XY( J) ='-XY ( 1 ) 

IF(CXX( J),GT„X)G0TP6P 
GOTO 61 
62 N=J-1 
KK = 1 

61 CONTI Nil F 
10? KK=0 
J = 0 
K = M 

00 71 1=1 , K 
IF(KK,EOol IGOTn?! 

IF(THX( I I, LT,,0 IGOTOTi 
J = J-H 

M = J 

THX( J)=THX( T ) 

THY( J» = THY( T ) 

IF(THX( I GTjXIGOTOTP 
G0TD71 

72 M=J-1 
^ — 1 

71 continue 

c. 

C determine frroR bar magmTTijoE 

c 

IF(N,fq, 0) GOTOl 03 
DO ■’00 1=1 

CTSDNi n = EXYm-nsopT(cxY(n''2co.i/RO"n 
CTsup( n = EXY( n+osoRT ( EXY( n '-P 00 . ) / ?on, 

700 CONTINOF 

c 

C INITIALIZE PLOT AND WRITE I DENT I f I CAT I OM 

c 

103 CALL PLOTS 

CALL symbol (0,'^, ,2°,TITLF0,C ,1 6 ) 

CALL PL0T(o,2„,-^» 

C 

C DRAW outline OE GRAPH 
C 

7 = 0 

CALL P|OT(7,Y,P) 

CALL PLOT(X,Y,P) 

CALL PLOT IX, 7,0) 

CALL PLOT! 7,7,0) 

C 

C DRAW OUTLINE OF TITLF BOX 
C 

Rl=-, ^ 

R7 = 4 ,4 

R3=?,S4 

B4=4.6 
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1 

H 



R5 = R! + > I 
B^=R3+, 1 

CALL t>L 0T( R1 , R^, 

CALL PLOT( q?, B3, 

CALL PI ny ( R? , R4, :> ) 

CALL P!.QT( R1 , R4, ->) 

CALL PI HT( fXl , RB,?) 

CALL SYMBOL ( R‘=; tR6, . 14, T1,0 

B5=R5+, 24 

CALL NnMRE'?(BB,R4,,14,T,0,0,3) 

RR = R R+1 „ 3? 

B6 = P6+-, 1 
cxP=7, 

CALL NDMRPP (R5,B6, •,07,cxP,0.'^,- 
RS^R'^+l , 70 
B6=B6- . 1 

CALL NUMBER (B5,B6, ^ 14, PT ,0,2 

B5 = 85-2,7/^ 

B6 — R^'^'', 24 

CALL SYMB0L(R5,R6, -,1.4, T?,0 

IF( N,EQ<.0) GOTOTOl 
R6=B6+o P? 

R5 = R5+-> B2 
R7 = B5+, R 

CALL PLOT (R5,R6,7) 

CALL PLOT (B7,RA,2) 

B6=B6-, 07 
B5 = B5--, R7 
701 R6=B6+^74 



CALL 


SYMBOL(R5,B6, ,14, 


T3,0 


B6=R6+, 24 




CALL 


SYMBOL (B5, R6,, 14, 


T4 ,0 



C PLP'T TIC MARKS 
C 

G=«! 

H= Y— o 1 

1570 call” PLOT(r,,Y, 3) 

CALL PL0T(G,H,2) 

G— G + « 1 

IP (GoLBnXlGO TO 1570 
G= X~ o 1 
H = ,l 

1580 CALL PLnT(G,H,3) 

CALL PL0T(G,O„,2) 

G= G~ « 1 

IF(gIgT«0, )G0 to 1 580 
G=l« 

H— 2 

1570 CALL PL0T(G,M,3) 

CALL PL0T(G,O,2) 

G=G+1, 

IF(G,LT,X) GO 70 1^70 
H=Y-,7 

1680 CALL PLnT(G,H,3) 

CALL PL0T(G,Y,7) 

G= G~ 1 o 

IF(’g,GT,0o) go to 16«0 
G = *l 
H — Y“ 1 

I'^Tl CALl’pLOT(G,H, 3) 

CALL PL0T(<'> ,H,2 ) 

H — H“ a 1 

TFIh’gT,^, ) GO TO 1571 

G=X-ol 

H= >I 

1501 CALL PL0T(G,H,3) 

CALL PL'77(X,H,2) 

H=H+,1 

1F(H«LT,Y) GO TO 1581 

G=X-,2 

H=Y-1, 



P,3?) 



1 ) 



0,37> 



0,32) 

0,32) 
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ooo non ooo ooo non ooo 



U>7] CALL PLnT(G,H,^) 

CALL PLnT(X,H,7) 

M” 1 <t 

IF(HoGT.O) to ]a-'1 

G=,P 

H=K 

16P1 CALL PLnT(G,H,A) 

CALL PLQT(0,H,P) 

H = H+U 

1F(HoLT,Y) go to 
LARFL AXIS 

CALL SYMROL (3,5,--, . 1 A , t iti c y , , 4 ) 

CALL SYWROL(-o?'=^,l,S6,,14,TTTlFY,oO,,?4» 

MO^^RFR energy axis 

1050 G=-,?P 
H = -« -^2 
FL0=SP 

l'A60 CALL NMMPPR(G,H,nl4,FLn,0,P) 

G=G+1> 

IF(G<,GT,X) GO TO 1070 
FL0=FLn + XT NCR 
GO TO 1060 

PLOT EXPERIMENTAL DATA 

1070 IF(N,EOoO)GOTniOA 

CALL LI NEC EXX , EXY, N, R , -6 ) 

PLOT ERROR flags 



DO 1010 I=1,N 
FS1 = FXX( n-,05 
*^S2=TXX ( I ) +.,05 
CALL PLOT! f^Sl ,CTSMP( I) ,R ) 
CALL PLOT( '"S?,CTSUD( I ) ,2 ) 
CALL PLO"^ ( FXX ( n ,CTSUP ( I ) ) 

CALL PL0T(EXX(n ,CTSnM( I ),R) 
call PLOT! ESI tCTSONI I ) ,0 ) 
CALL PLOT! '=SR,CTSDM( I) ,2) 

1010 continue 

PLOT THEORY CURVE 

104 CALL LINE! THX,THY,m, 2, n 

FINALIZE PLOT 

CALL PLOTIO. ,P,3 ,-3) 

CALL PLOTF 

return 

END 
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